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Abstract 



The degree of isospin fractionation is measured by (N/Z) n / (N/Z) N . mf , where 
(N/Z) n and (N/Z) ^ imf are the saturated neutron- proton ratio of nucleon emissions ( 
gas phase) and that of fragment emissions (liquid phase) in heavy ion collision at in- 
termediate energy . The calculated results by using the isospin-dependent quantum 
molecular dynamics model show that the degree of isospin fractionation is sensi- 
tive to the neutron-proton ratio of colliding system but insensitive to the difference 
between the neutron-proton ratio of target and that of projectile. In particular, 
the degree of isospin fractionation sensitively depends on the symmetry potential. 
However its dependences on the isospin dependent in-medium nucleon-nucleon cross 
section and momentum dependent interaction are rather weak. The nucleon emis- 
sion (gas phase) mainly determines the dynamical behavior of the degree of isospin 
fractionation in the isospin fractionation process, compared to the effect of frag- 
ment emission. In this case, we propose that (N/Z) n / (N/Z) N . mf or (N/Z) n can be 
directly compared with the experimental data so that the information about sym- 
metry potential can be obtained. 
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1 Indroduction 

The studies of the isospin effects at intermediate energy heavy ion collision can be 
used to get the information about isospin dependent in-medium nucleon-nucleon 
cross section and isospin dependent mean field (symmetry potential). In order 
to obtain this information several interesting isospin effects in heavy ion colli- 
sions have been explored both experimentally and theoretically over the last few 
years[I — 16,26 — 28]. Bao-An Li et al investigated the isospin effect of the mean 
field and showed that the neutron-proton ratio of preequilibrium nucleon emis- 
sions, the neutron-proton differential flow and proton ellipse flow are probes for 
extracting the information of the symmetry potential [1,2,14,15,16]. R.Pak et 
al have found that the isospin dependences of collective flow and balance energy 
are mainly originated from the isospin-dependent in-medium nucleon-nucleon cross 
section [17,18]. In order to extract above infoirmation Bao-An Li ,M.Di Toro ,M. 
Colonna and V.Baran et al studied the isospin fractionation in the intermediate 
energy heavy ion collisions in recent years [15, 16, 19,30]. The isospin fractionation 
is an unequal partitioning of the neutron to proton ratio N/Z of unstable isospin 
asymmetry nuclear matter in the low and high density regions. Bao An Li 's work 
proposed a very useful point of review about isospin fractionation [14]. V.Baran et 
al studied the isospin fractionation dynamics in details and gave more interesting 
results to explain the isopin fractionation dynamics [30]. An indication of this phe- 
nomenon has been found recently in the intermediate energy heavy-ion experiments 
by H. Xu et al [9]. 

However the two essential ingredients in heavy ion collision dynamics, the isospin de- 
pendent in-medium nucleon-nucleon cross section and the symmetry potential have 
not been well determined so far. In order to compare with experimental data directly 
it is important to quantify the gas phase and the liquid phase as well as their depen- 
dences on the dynamical ingredients in isospin fractionation process. Based on the 
isospin dependent quantum molecular dynamics model (IQMD)[2, 20, 21, 22, 26] we 



investigated the degree of isospin fractionation (N/Z) n /(N/Z)N imf , where (N/Z) n 
and (N/Z)Ni m f are the saturated neutron to proton ratio of nucleon emissions and 
that of the fragment emissions respectively. We investigated the dependeces of 
(N/Z) n /(N/Z) Nim f on the symmetry potential, isospin dependent in-medium nu- 
cleon nucleon cross section and momentum dependent interaction in some reaction 
processes. The calculated results show that (N/Z) n /(N/Z)mmf only sensitively de- 
pends on symmetry potential for neutron-rich colliding system. We simulated and 
discussed the dependences of isospin fractionation degree on the neutron proton ratio 
of colliding system and the difference between neutron proton ratio of projectile and 
that of target. The roles of gas phase and liquid phase in the isospin fractionation 
process are also investigated in detail. 

2 IQMD Model 

As we have known that quantum molecular dynamics (QMD)[21, 22] contains two 
dynamical ingredients: density dependent mean field and in-medium nucleon nu- 
cleon cross section. To describe isospin effects appropriately, QMD should be mod- 
ified accordingly. The density dependent mean field should contain correct isospin 
terms, including the symmetry potential and the Coulomb potential. The in-medium 
N-N cross section should identify neutron-neutron , proton-proton and neutron- 
proton collisions, in which Pauli blocking should be counted by distinguishing neu- 
trons and protons. 

Considering the above ingredients, we made important modifications in QMD to ob- 
tain the isospin dependent quantum molecular dynamics (IQMD) . In IQMD model 
the initial density distributions of the colliding nuclei are from the calculations by 
applying the Skyrme-Hatree-Fock model with the parameter set SKM*^. The 
code of IQMD without the collision term is used to obtain the ground state proper- 
ties of the colliding nuclei. The ground state properties, such as binding energies and 
RMS radii etc, are consistent with the experimental data. Thus the parameters of 
interaction potentials are fixed and used as input data in the dynamical calculations 
by IQMD. In our calculations the interaction potentials are taken as: 

U(p) = u Sky + jj Coul + JJ sym + JJ Yuk + JJ MDI + jj Pauli (1) 



where U Sky , JJ Coul , JJ Yuk and JJ Pauh are Skyrme potential, Coulomb potential,Yukawa 
potential and Pauli potential. JJ MDI is momentum dependent interaction with the 
form of 

U MDI = tJn 2 Mpl~ti) 2 + l]-- (2) 

Po 

The more detailed physics ingredients and their numerical realization in the IQMD 
model can be found in Refs [Chapter 10 in Ref.2,20,21,22,26]. U sym is the symmetry 
potential. In this paper, three different forms of JJ sym are used [1, 15], 

Ur = cF 1 (u)5r z (3) 

U° ym = cF 2 (u)[5r z - l -5 2 } (4) 
_ 2cF 3 {u)5t z cF 3 (u)S 2 

Us ~ i + u + ~(TTW {) 

with 

1 1 for neutron 



z 



— 1 for proton 

Here c is the strength of symmetry potential, with the value of 32MeV, Fi(u)=u 
, Fo(u) = u 1 / 2 and FJu)=u 2 , u = —. 5 is the relative neutron excess 5 = Pn ~ Pp = 
Pn ~ pp . p ) p o) p n and p p are the total , normal , neutron and proton densities respec- 
tively. 

The curvatures K sym of U{ ym , U s 2 ym and U^ ym are -27 MeV, -69 MeV and 61 MeV 
respectively. 

An empirical density dependent expression of the in-medium N-N cross section is 
[24]: 

cnn = (1 + a—)crfj^ e (6) 
Po 

with the parameter a ~ —0.2 which can well reproduce the flow data . Here crjj^ is 
the experimental nucleon nucleon cross section[25]. The free neutron-proton cross 
section is about 3 times larger than the free proton-proton or the free neutron- 
neutron cross section below about 400 MeV. It should be pointed out that the ratio 
tfnp/o'pp in the nuclear medium sensitively depends on the evolutions of the nuclear 
density distribution and beam energy. Here we made use of equation (6) to take the 
medium effect of two-body collision into account, in which the neutron-proton cross 
section is always larger than the neutron-neutron or proton-proton cross section in 



the medium at the beam energies studied in this paper. 



3 Results and Discussions 

The isospin effect of the in-medium nucleon nucleon cross section on the observ- 
ables is defined by the difference between the observables from an isospin dependent 
nucleon nucleon cross section a lso and those from an isospin independent nucleon 
nucleon cross section a nmso . The a lso is illustrated by a np > cr nn —a pp and a nmso 
means a np = a nn = a pp where a np , <J nn , a pp are neutron proton, neutron neutron and 
proton proton cross sections in medium respectively. In (N/Z) n /(N/Z) N . mf , the 
multiplicity of intermediate mass fragments N im f includes all fragments with the 
charge number from 2 to (Z p + Z t )/2 , where Z p and Z t are the charge numbers of 
the projectile and target respectively . 



3.1 The dependences of the degree of isospin fractionation 
on the symmetry potential and isospin dependent im- 
medium nucleon nucleon cross section 

Fig.l shows the impact parameter averaged values of < (N/Z) n > b / < (N / Z) N . mf>b 
as a function of the beam energy E in the reaction system 124 Sn + 124 Sn. There 
are four different cases: JJ{ ym + a lso (solid line), U^" 1 + cr tso (dashed line), £/| ym + 
<r noiso (dot-dashed line) and U{ ym + a noiso (dot line). 

It is found that for the different nucleon nucleon cross sections, a lso and a nmso , 
the variations among the (N/Z) n /(N/Z) im f are smaller when the symmetry poten- 
tial is fixed as either C/| ym or Ul ym . However for different symmetry potentials U^ m 
and JJ{ ym the gaps between the (N/Z) n /(N/Z) im f are larger when the nucleon nu- 
cleon cross section is fixed as either a lso or a nmso . So it might be concluded that 
(N/Z) n /(N/Z) im f sensitively depends on the symmetry potential and weakly on 
the isospin effect of in-medium nucleon nucleon cross section. More complete un- 
derstanding of the isospin fractionation processes can be obtained from the analysis 
of the density dependent chemical potentials for neutrons and protons [30]. As one 



part of chemical potential, the symmetry potential is repulsive for neutrons and at- 
tractive for protons in the expanding process of the colliding system. This leads to 
the different chemical potential degradations with nuclear density p q for the neutron 
and proton in the low density region below normal nuclear density. However the 
mass flow in the nonequilibrium expanding process is determined by the difference 
in the local values of the chemical potential and such difference directs the system 
from higher chemical potential region into lower values until equilibrium. This trans- 
port effect induces the unstable isospin asymmeric nuclear matter to be separated 
into a neutron-rich low density phase and a neutron-poor high density phase for the 
neutron-rich system. Namely it is the different chemical potential degradations with 
nuclear density for neutrons and protons below normal nuclear density that induces 
the isospin fractionation. 

According to above analysis the symmetry potential in the chemical potential is a 
main dynamical ingredient for producing the isospin fractionation but the influence 
of two-body collision on it is rather weak in the beam energy region in this paper. 

3.2 The dependence of [N jZ) n j{N jZ)^ imf on neutron to pro- 
ton ratio of the colliding system 

In order to investigate the evolution of (N/Z) n /(N/Z) N . mf with the neutron to 
proton ratio of the colliding system , Fig. 2 shows the (N/Z) n /(N/Z) N . mf as a func- 
tion of the neutron-proton ratio of the colliding system for the reactions 104 Sn + 104 
Sn, n2 Sn + 112 Sn and 124 Sn + 124 Sn at the beam energy E = 50 MeV/nucleon , 
impact parameter b=0.0 fm (right window) and 4.0 fm (left window) . The neu- 
tron to proton ratios of the above three colliding systems are 1.01,1.34 and 1.48. 
In Fig. 2 there are three different symmetry potentials U{ ym , U^" 1 and U^ ym . The 
calculation results of the three reactions show that the (N / Z) n / (N / Z) x imf and its 
dependence on the symmetry potential are enhanced with increasing the neutron to 
proton ratio of the colliding systems. In particular, the (N/Z) n /(N/Z) N . mf is larger 
than 1.0 for the neutron-rich system 124 Sn + 124 Sn , which means the neutron-rich 
gas phase and the neutron-poor liquid phase , compared to the neutron to proton 
ratio of colliding system. It is also found that the [N / Z) n / {N / Z)^. mf is enhanced 



with increasing symmetry potential strength. The amplitudes of (N/Z) n / (N/Z)N imf 
from large to small correspond to the order of symmetry potential strength. The 
order of symmetry potential strengthes in the lower density region below the normal 
nuclear density are from U s 2 ym (dot line) , U{ ym (solid line) to U^ ym (dot-dashed line) 
for the neutron-rich system. The order of symmetry potential strength below the 
normal density is contrary to those at the normal nuclear density (see Fig. 4.1 in [1]). 
However for the neutron-poor systems, such as, 104 S , n+ 104 Sn and 112 Sn+ 112 Sn the 
values of (N/Z) n /(N/Z) N . mf are decreased gradually below 1.0 with decreasing the 
neutron-proton ratio of colliding system. This means the neutron-poor gas phase 
and neutron-rich liquid phase, compared to the neutron proton ratio of colliding 
system. In the isospin fractionation process for the neutron-poor system, (N/Z) n 
is decearsed due to neutron-poor colliding system. On the other hand , (N/Z) im f 
is increased, compared to the initial neutron proton ratio of neutron-poor colliding 
system. This can be explained by the fact that the excited primary neutron-poor 
fragments decays to a stable neutron-rich one which is closer to the beta stable line. 

3.3 The dependences of < (N/Z) n > b / < (N/Z) N . mf > b on the 
gas phase and liquid phase 

In order to study the dependences of < (N/Z) n > b / < (N/Z) N . mf > b on the 
gas phase and liquid phase, Fig. 3 shows the impact parameter averaged values of 
< (N/Z) n > b (in the left panel) and < (N/Z) N . mf > b (in the right panel) as a 
function of the beam energy E for the different symmetry potentials U{ ym and U^™- 
The incident channel conditions and the line symbols in Fig. 3 are the same as those 
in Fig.l. It is clear to see that the dynamical behavior of < (N/Z) n > b (gas phase) 
is very similar to that of < (N / Z) n > b / < (N/Z) Nimf > b in Fig.l. However 
all lines for < (N/Z)N imf > b (liquid phase) corresponding the different symmetry 
potentials JJl ym and U^ ym as well as different nucleon nucleon cross sections a lso 
and a nmso are very close to each other. Namely < (N/Z)N imf > b is insensitive to 
both the symmetry potential and the isospin effect of in-medium nucleon nucleon 
cross section. So the nucleon emissions (gas phase) mainly determine the dynamical 
behavior of the degree of isospin fractionation for the neutron-rich colliding system 
in the beam energy region studied here. 



It should be pointed out that the calculated results in above three figures are taken 
by accounting the time up to 300 fm/c. This time roughly corresponds to the freeze- 
out time because during fragmentation stage the colliding system breaks up and the 
fragment formation process occurs up to the time around 240-300 fm/c in the beam 
energy region studied here. From the time evolusion of impact parameter averaged 
value of nucleon emission number < N n > b and that of fragment emission number 
< N im f >b in the reaction 124 Sn + 124 Sn at E= 50 MeV/nucleon shown in Fig.4.b, 
we can see that < N im f > b almost stabilizes after 250 fm/c, even though < N n > b 
continues to be enhanced very slowly after 250 fm/c. However the impact parameter 
averaged value of the neutron proton ratio of nucleon emission and that of fragment 
emission < (N/Z) n > b and < {N/Z)^ imf > b shown in Fig. 4. a are almost constant 
after 200 fm/c. 

From above analysis and discussion in Fig.4 we might conclude that even though 
the simulation by IQMD could not include all secondary decaies, it does not matter 
too much on the values of < (N/Z) n > b and < (N/Z) N . mf > b at 300 fm/c and the 
conclusions. Furthermore the observable for comparing with the experimental data 
in this paper is < (N/Z) n > b /< (N/Z) N . mf > b . 

3.4 The dependence of (N/Z) n /(N/Z)]y. mf on the difference 
between neutron to proton ratios for projectile and tar- 
get 

Fig.5 shows the time evolution of (N/Z) n /(N/Z) N . mf in the reaction systems which 
have different neutron-proton ratios for target and projectile but with the same mass 
and the same neutron-proton ratio of colliding system ,the same JJl ym and the same 
a %so at beam energy E=50 MeV/nucleon and impact parameter of 2.0 fm . In Fig.5 
there are two couples of reactions 74 Zn + 74 Se and 74 Ge + 74 Ge ( left window) as 
well as 132 Sn+ 112 Sn and 124 Sn+ 124 Sn (right window). The neutron-proton ratios 
of projectile, target and colliding system are listed in table 1. 
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From Fig. 5 it is clear to see that all solid lines and all dashed lines are very close 
to each other for the same mass and the same neutron proton ratio of each couple 
of colliding systems. It shows that the difference between the neutron-proton ratios 
for target and projectile doesn't influence isospin fractionation process, instead it 
sensitively depends on the neutron-proton ratio of colliding dsystem as observed in 
Fig.2 

It is worth mentioning that because the scaling time of the isospin fractionation 
process is faster than that of fragmentation process , (N/Z) n is much larger than 
(N/Z) N . mf before 125 fm/c . This indicates that the large fluctuation before 100 
fm/c, then (N/Z) n and (iV/Z)jv im/ gradually approach their equilibrium values with 
the increasing colliding time. This is why the beginning time of the time evolutions 
for figures is started at 100 fm/c. 



3.5 The role of MDI on the dynamical process of isospin 
fractionation 

The nonlocal property of the nuclear interaction leads to a repulsive momentum 
dependent interaction (JJ MDI ) in the intermediate energy heavy ion collisions. In 
order to study the influence of the JJ MDI on the dynamical process of isospin frac- 
tionation Fig. 6 shows the time evolution of impact parameter averaged values of 
< (N/Z) n > b I < (N/Z) N . mf > b with the same a lso but different symmetry poten- 
tials Ul ym and U^™ ' . In calculations we consider momentum dependent interaction 
U MDI and momentum independent interaction jj noMDI for the reaction 12A Sn+ 124 Sn 
at E= 50 MeV/nucleon. There are four cases: U[ ym +U MDI ; U° ym +U noMDI ;C/f m 
+U MDI and U s 2 ym + U noMDI . The gap between the values of < (N/Z) n > b / < 
(N/Z) N . mf > b is smaller for the U MDI and U noMDI with the same U° ym or XJ sym . 
However, the gap between the values of < (N/Z) n > b / < [N / Z)^ irnf > b is lager for 
the same the JJ MDI or jj noMDI but different symmetry potentials. It turns out that 
the influence of JJ MDI on the isospin fractionation process is not obvious because 
the emission probabilities induced by JJ MDI are about the same for the neutrons 
and protons . 



4 Summary and conclusions 

In summary, from the calculation results we can get following conclusions: (l)The 
< (N/Z) n >b / < (N/Z)N imf >b is enhanced with increasing the neutron to proton 
ratio of colliding system and symmetry potential strength for neutron-rich system. 
However, the isospin fractionation process is insensitive to all these ingredients for 
the neutron-poor system in beam energy region studied here. 

(2) In particular^ (N/Z) n > b / < (N/Z) N . mf > b and < (N/Z) n > b sensitively 
depend on the symmetry potential but weakly on the isospin effect of in-medium 
nucleon-nucleon cross section and momentum dependent interaction for the neutron- 
rich system. 

(3) The dynamical behavior of < (N/Z) n > b / < (N/Z) N . mf > b is mainly determined 
by the nucleon emissions (gas phase). 

(4) The (N/Z) n /(N/Z)N imf sensitively depends on the neutron to proton ratio of 
colliding system but weakly on the difference between the neutron proton ratios for 
target and projectile. 
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Figure captions 

Fig.l The impact parameter averaged value of < (N/Z) n / < (N/Z)^ imf>b as a 
function of the beam energy E for the reaction system 124 Sn + 124 Sn and the 
different symmetry potentials U{ ym and XJ s 2 ym with a iso or a noiso . 

Fig. 2 The (N/Z) n / [N/Z)^ im} as a function of the neutron-proton ratio of colliding 
systems 104 Sn + 104 SnJ l2 Sn + 112 Sn and 124 Sn + 124 Sn for three different sym- 
metry potentials U{ ym ,XJ^ m and U^ ym at the beam energy of 50 MeV/nucleon 
and impact parameter of 4.0 fm (left panel) and 0.0 fm (right panel). 

Fig. 3 The impact parameter averaged values of < (N/Z) n > b (left panel) and 
< (N/Z) N . mf > b (right panel ) as a function of beam energy E for the different 
symmetry potentials U{ ym and C/| ym as well as the reaction 124 Sn + 124 Sn. 



Fig. 4 The time evolusions of impact parameter averaged values of nucleoli emission 
number < N n > b and fragment emission number < N im f > b (Fig.4.b)) as well 
as that of neutron proton ratio of nucleon emission < (N/Z) n > b and that of 
fragment emission < (N/Z) N . mf > b (Fig. 4. a) for the reaction 124 Sri + 124 Sn at 
E= 50 MeV/nucleon. 

Fig. 5 Time evolution of (N/Z) n /(N/Z) N . mf for different neutron-proton ratios for 
target and projectile but the same mass and the same neutron proton ratio of 
each couple of colliding systems for the reactions 74 Zn + 74 Se and 74 Ge + 74 Ge 
(left window) as well as 132 Sn + 112 Sn and 124 Sn + 124 Sn (right window). 

Fig. 6 the time evolutions of the impact parameter averaged values of < (N/Z) n > b 
I < (N/Z) N . mf > b with the same a iso ,the same U MDI or U noMDI but different 
symmetry potentials JJl ym and U^™ for the reaction 124 Sn+ 12A Sn at the beam 
energy of 50 MeV/nucleon. 



